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ABSTRACT 

We show that if in a classical accretion disk the thin disk approximation fails interior 
to a certain radius, a transition from Keplerian to radial infalling trajectories should 
occur. We show that this transition is actually expected to occur interior to a cer- 
tain critical radius, provided surface density profiles are steeper than cx 
and further, that it probably corresponds to the observationally inferred phenomena 
of thick hot walls internally limiting the extent of many stellar accretion disks. Once 
shears stop, the inner region of radially infalling orbits is naturally expected to be cold. 
This leads to the divergent focusing and concentration of matter towards the very cen- 
tral regions, most of which will simply be swallowed by the central star. However, if a 
warm minority component is present, we show through a perturbative hydrodynam- 
ical analysis, that this will naturally develop into an extremely well coUimated and 
very fast moving pair of polar jets. A first analytic treatment of the problem described 
is given, proving the feasibility of purely hydrodynamical mechanisms of astrophysi- 
cal jet generation. The purely hydrodynamic jet generation mechanism explored here 
complements existing ideas focusing on the role of magnetic fields, which probably 
account for the large-scale coUimation and stability of jets. 
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1 INTRODUCTION 

Astrophysical jets occur over a large range of astrophysi- 
cal scales, from the stellar Newtonian scales of HH objects 
(e.g. Reipurth & Bally 2001), to the relativistic cases of 
microquasars and gamma ray bursts (e.g. Mirabel & Ro- 
driguez 1999, Meszaros 2002), to the Mega-parsec extra- 
galactic scales of AGN jets (e.g. Marscher et al. 2002). Al- 
though the processes of propagation and coUimation appear 
to be relatively well understood in terms of the interplay 
between the hydrodynamics of the problem (Scheuer 1997) 
and magnetic fields (e.g. Blandford 1990), the precise mech- 
anism of jet generation probably remains as the most uncer- 
tain part of the system. 

Existing jet generation mechanisms focus on the inter- 
action between the rotating material in the inner regions 
of the accretion disk and the magnetic field of the central 
(rotating or otherwise) star, or perhaps an external mag- 
netic field threading the disk. See for example Blandford & 
Payne (1982), Henriksen & Valls-Gabaud (1994), Meier et 
al. (2001), Price et al. (2003) for some examples. Although 



generically successful, detailed comparisons with observa- 
tions do dot always yield consistent results, e.g. Ferreira et 
al. (2006). Further, the details of relative orientation of stel- 
lar spin, accretion disk and magnetic field configurations, 
have to be supplied in somewhat highly specific manners. 

As an alternative, we explore the possibility of hydro- 
dynamical jet generation mechanism where magnetic fields 
play no role, but only the intrinsic hydrodynamical physics 
of the interplay between the accretion and the central star 
wholly determine the characteristics of the jet. The main 
obstacle to such a scheme is the presence of a centrifugal 
barrier associated with the angular momentum content of 
the material in the accretion disk. We show however, that 
the thin disk approximation for accretion disks, which is 
equivalent to the assumption of quasi-circular orbits for the 
disk (e.g. Pringle 1981), is expected to break down inter- 
nal to a certain critical radius, resulting in a transition to 
quasi-radial flow for the disk. Once this obstacle has been re- 
moved, an analytic first order perturbation treatment serves 
to demonstrate the feasibility of purely hydrodynamical jets. 
We show also that many generic features of observed astro- 



2 X. Hernandez, P. L. Rendon, R. G. Rodriguez-Mota, A. Capella 



physical jets across a wide range of scales, can be naturally 
accounted for in the general model presented. The presence 
of magnetic fields in jet phenomena is evident empirically, 
however, it is not impossible that their main contribution 
to the problem could be restricted to the radiation of the 
jet material and to its long range coUimation, with purely 
hydrodynamical physics playing a part in the actual jet gen- 
eration mechanisms. 

The problem of jet formation has been studied exten- 
sively in the context of classical hydrodynamics, most often 
regarding fluid-body interactions. The appearance of stable 
coaxial jets resulting from radially-symmetric velocity fields 
over thin fluid sheets has been established by, among oth- 
ers, Taylor (1960). The role played by both the magnitude 
and the direction of velocity in the formation of this type of 
jet is the subject of a theoretical study by Glauert (1956), 
where it is shown that at the point at which the jet forms, 
a large velocity gradient is observed, and momentum flux 
is constant, with horizontal momentum being transformed 
into vertical momentum. Similar results were obtaiucd by 
King & Needham (1994), who provide an asymptotic study 
of a jet formed by a vertical plate accelerating into a semi- 
infinite expanse of stationary fluid of finite depth with a free 
surface and a gravitational restoring force. It is found that 
as the fiuid approaches the plate in a horizontal direction, 
a gradual rise in free-surface elevation occurs. Eventually, a 
thin region is reached where vertical velocity dominates hor- 
izontal velocity as a consequence of the fiuid finding it more 
difficult to overcome the inertia it would meet by continuing 
its horizontal motion than escaping vertically towards the 
low-pressure free surface. In essence, this same mechanism 
can allow for jet formation in the context of the problem we 
study hero. 

In section 2 we develop the criterion for transition to 
radial flow in a standard accretion disk. The perturbative 
solution to the resulting problem of a radially infalling disk 
is then developed in section 3. Section 4 presents trajectories 
for particular cases of the solution obtained in the previous 
section, in dimensionless units. Finally, we give our conclu- 
sions in section 5, 



2 THE TRANSITION TO RADIAL FLOW 

We start from a standard thin accretion disk where material 

orbits on quasi-circular Keplerian orbits around a central 
star of mass M. Assuming axial symmetry and cylindrical 
coordinates, the total angular momentum of a ring at radius 
R of width A_R wiU be: 



L = 27vR-'ARE{R)n{R) 



(1) 



where E and Q are the surface density and orbital frequency 
profiles of the disk, respectively. If the breaking torque on a 
given ring of the disk, due to its exterior ring is denoted by 
t{R), the total breaking torque on the ring at radius R will 
be: 



TB =t{R + AR)-t{R). 



(2) 



The rate of change of the angular momentum of the ring 
in question can now be calculated as: 



tb 



(3) 



If at any radius a substantial fraction of the angular 
momentum of the ring is lost through viscous torques over 
an orbital period, the assumption of quasi-circular orbits will 
break down, and the disk will make a transition to mostly 
radial orbits. This condition can be stated as: 

^ < n • 

Substitution of equations (1) and (2) into the above 
yields: 



R^Q^T. < 2 



(5) 



as the condition for the onset of radial flow in the disk. 

We can now introduce a model for the torques in terms 
of the rate of shear and the effective viscosity ly (e.g. Von 
Wiezsacker 1943, Pringle 1981) as: 

'dR' 

Taking = GM/R^ to substitute the above equation 

into condition (5) gives 



(6) 



dR 



(7) 



To proceed further we can take for example, 
v =constant and a model for the disk surface density profile 
of the form: 



(8) 



of the type often used in models of accretion disks when 
fitted to observations (e.g. Hughes et al. 2009). Use of the 
two above forms for v and E(ii) reduces condition(7) to: 



6w(n- 1/2) > {GMRf^. 



(9) 



The above condition will always be met in any accretion 
disk with n > 1/2, interior to a transition radius Rt given 
by: 



Rt — 



[67rt/(n- 1/2)]^/^ 



GM 



(10) 



It is interesting that directly observed accretion disks 
have spectra which when modelled typically yield n ~ 1, 
in general < n < 1.5, e.g. Hartmann ct al. (1998), Lada 
(2006), Hughes et al. (2009). This leads one to expect the 
transition to radial fiow to occur in many of the stellar ac- 
cretion disks, interior to radii given by equation (10). 

Finally, we can explore the consequences of introducing 
an a prescription (Shakura & Sunyacv 1973) into condition 
(9), V = acH, where c and H are respectively the sound 
speed and height of the disk, with a a dimensionless number, 
yielding the dimensionless condition: 

Al(J) <677a(n-l/2), (11) 

where M is the ratio of the Keplerian orbital velocity in the 
disk to the sound speed in the disk. Although the above con- 
dition is only valid for constant u, it illustrates the equiv- 
alence between the assumption of quasi-circular orbits for 
the material in the disk, and the assumption that the disk 
is thin. We see that the breakdown of the assumption of 
quasi-circular orbits (condition 9), corresponds to the point 
where disk is fat. 

Writing M in astrophysical units as: 
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V T 



M 



0.5Mfr 



1/2 



(12) 



we see that for values typical of what is observationally in- 
ferred for the stellar mass and position and temperature of 
accretion disk walls in T Tauri protoplanetary disks, those 
indicated in the above equation (e.g. D'Alessio et al. 2005, 
Espaillat et al. 2007, Hughes et al 2009), one should expect 
the breakdown of the thin disk approximation {M » 1, 
e.g. Pringle 1981), and consequently, the transition to radial 
flow. 

Taking typical inferred values for R/H at the wall of 
~ 5 (e.g. D'Alessio et al. 2005, Espaillat et al. 2007, Hughes 
et al 2009) , we can now write the condition for the transition 
to radial flow, locally at the wall, as: 



10 = 67ra(n-l/2). 



(13) 



We see that for a standard value of n ~ 1 the above 
equation implies a reasonable value of a ~ 1 at the thick 
wall, significantly higher than the values of ^ 0.01 and lower, 
which apply for the body of the thin accretion disk beyond 
this radius. A substantial increase of a as R/H decreases is 
expected in any turbulence driven viscosity model for accre- 
tion disks, e.g. Firmani, Hernandez & Gallagher (1996), in 
the context of galactic disks. 

In terms of the debate surrounding the inference of inner 
holes in observed accretion disks, many solutions have been 
proposed in terms of disk clearing mechanisms; grain growth 
(e.g. Strom et al. 1989, Dullemond & Dominik 2005), pho- 
toevaporation (e.g. Clarke et al. 2001), magnetorotational 
instability inside-out clearing (e.g. Chiang & Murray-Clay 
2007, Dutrey et al. 2008), binarity (e.g. Ireland & Kraus 
2008) and planet-disk interactions (e.g. Rice et al. 2003). 
None of the above is entirely satisfactory, as noted by Hughes 
et al. (2009), mostly due to their incompatibility with a 
steady state solution. An alternative solution under the pro- 
posed scenario, is that there is no actual disk material clear- 
ing, only a transition to radial flow at the thick wall, and 
consequently a shear- free flow interior to this point. Once 
the disk heating mechanism is removed, as one should ex- 
pect from the analysis presented in this section, the inner 
disk disappears from sight. 



3 HYDRODYNAMICAL JET SOLUTIONS 

We shall now model the physical situation resulting from 
the scenario described above as an axially symmetrical dis- 
tribution of cold matter in free fall towards a central star of 
mass M. Taking a spherical coordinate system with 9 the 
angle between the positive vertical direction and the position 
vector rwe have: 



1 dir'^pV) 

, u dv 



d{sin{e)pU) 



rsin{6) 



0, 
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1 dP _ GM 
1 dP 

Vp~de' 



(14) 
(15) 
(16) 



for the continuity equation, and the radial and angular com- 
ponents of Euler's equation. In the above, V , U, p and P 



give the radial velocity, angular velocity, matter density and 
pressure, respectively. We have neglected temporal deriva- 
tives, as we are interested at this point, in the characteristics 
of steady state solutions. We take as a background state a 
free-falling axially symmetrical distribution of cold gas, de- 
scribed by Vo = -(2GM/r)^/^ Uo = and VPo = 0, a 
consistent solution to eqs.(15) and (16). Having ignored the 
inclusion of the radius at which the transition to radial flow 
takes place in the choice of Vo limits the validity of the anal- 
ysis to radial scales along the plane of the disk rrmch smaller 
than Rt- This is justified by the fact that jets appear as phe- 
nomena extremely localised towards R 0. We now take a 
density profile given by: 



po{r,e) = f{r)gie), 



(17) 



were /(r) is a dimensionless function of r and g{6) describes 
the polar angle dependence of the infalling material, for ex- 
ample, one can ask for g{9 — it/2) — po, diminishing sym- 
metrically towards the poles. The choice of this last function 
will determine the details of the problem, and can be thought 
of as something of the type 



g{e) = poe 



\ V2eo J 



(18) 



with po a normalisation constant and 6o a form constant 
describing the flattening of the disk of infalling material. 
However, we shall mostly leave results indicated in terms of 
g{9). The continuity equation (14) now fixes /(r) through: 

g{9) d 



dr 



[(2GM)'''V^/VW] = 0, 



(19) 



and hence f{r)r^^^ = de., which completes the description 
of the background state through: 



Po(r,e) = ( 



9(0), 



(20) 



with r a constant which determines the point at which g{9 = 

7r/2) ^ Po = po- 

We shall now analyse the behaviour of a fraction of 
somewhat hotter material, through a first order perturba- 
tion analysis of the above solution. This small fraction could 
represent the last of the material to fully cool, or result from 
the irradiation of the central star onto the upper and lower 
surfaces of the radially infalling disk, as described by e.g. 
HoUenbach et al. (1994) or Alexander et al. (2006) for the 
standard case of Keplerian accretion disks, in connection 
with the problem of disk photoevaporation. In the above it 
is shown that ionizing radiation from the star generically 
creates an ionized layer on the surface of the disk. 

Regarding the evolution of this component, we shall also 
be interested interested in a steady state solution given by: 



Vir,9) = Vo{r) + Vi{r,9), 
U{r,9)^0 + Ui{r,9), 
p{r,9)=po{r,9)+pi{r,9), 



(21) 
(22) 
(23) 



where quantities with subscript (1) denote the perturbation 
on the background solution. Writing eqs.(14), (15) and (16) 
to first order in the perturbation one obtains after rearrang- 
ing terms: 

3(e)^lli^!^-B? - i)(^u,(0)!j{0)[A) 



dr 



dr 



r^/'^sin{9) d9 



4 X. Hernandez, P. L. Rendon, R. G. Rodriguez-Mota, A. Capella 



d {Vi/r 



l/2\ 



dr 



Ar'l'^ dpi 
g{0) dr- 



d{rUi) _ Ar' dpi 



dr 



gio) de ' 



(25) 



(26) 



In the above we have assumed an isothermal equation of 
state for the perturbation Pi = c^pi, an assumption often 
used in the modeling of astrophysical jets, e.g. the T Tauri 
jets observed and modelled by Hartigan et al. (2004). This 
idealised case serves to illustrate clearly the consequences of 
the physical setup being considered, as it allows for an ana- 
lytic solution. The generalisation to more realistic adiabatic, 
polytropic or otherwise equations of state can be performed 
numerically, and can be expected to yield qualitatively simi- 
lar results, although interesting differences in the details can 
be expected to emerge, which will be considered latter. In 
the above three equations we have introduced the constants 
A = (cV2GMf3)i/2 and B = {2GM/f''Y^^. 

To make further progress we can attempt a solution 
through the method of separation of variables, proposing a 
solution of the form: Vi = VrVe, Ui — UrUg, pi — prpe- This 
ansatz yields two independent systems of three equations 
each, one for the radial, and one for the angular dependences 
of the perturbations. The radial equations become: 



dr 

djrUr) 
dr 



dr 



Cgr-^pr 



d{r'/^Vr) /BCr 
dr \ A 



d{r' 



3/2 



pr) 



dr 



Cc 



Ur 

ri/2 ■ 



While the angular ones result in: 

A 



Pe, 



l{sin{e)g{e)Uo) 

de 



= -Ccsin{e)g{e)V0. 



(27) 

(28) 

(29) 

(30) 
(31) 

(32) 



In splitting the radial and angular dependences of the 
perturbed continuity equation, eq.(24), we have used the 
result of the angular equation of the perturbed radial Eu- 
ler equation, eq.(30). The constants Cr, Cg and Cc are the 
separation constants of the problem. Firstly we turn to the 
angular system, which allows an exact solution. 

We can take eq.(32) and substitute into it the product 
g{9)Ue from eq.(31), and the product g{0)Ve from eq.(30) to 
obtain an equation involving only pg: 



having solution: 



fCcCg\ ^ 



(33) 



pe = ciPm{cos{9)) + C2Q,n{cos{0)). 



In the above equation ci and C2 arc normalisation con- 
stants, and Pm and Qm are the Lcgcndrc polynomials of the 
first and second kinds, respectively. The index of these func- 
tions is given by the relation 2m = (4Cr -|- 1)^''^ — 1, where 
Ct = (CcCe/Cr). As typical of separation of variables prob- 
lems, we see the solution as an eigenvalue problem, with the 



separation constants determining the order of the solution 
function. The requirement of axial symmetry forces m to 
be even. Any such desired angular distribution for pg can 
now be constructed as an infinite series of the above func- 
tions. For simplicity we analyse the case of m = 2 (Ct = 6), 

Ci = pg, C2 = 0: 



pg = pg{3cos {6) - 1) 



(35) 



This will result in slightly less material along the plane of 
the disk, and slightly more along the poles, for the total 
density p = po + pi, with respect to the background state 
po- At this point eqs.(30) and (31) can be used to obtain: 



Vg=(^'^ /3cO£^(e)_l 



Ug 



f &Apg \ 
V Cg J 



g{0) 



cos(9)sin{9) 



g{0) 



(36) 



(37) 



The case of the radial system is more cumbersome, as 
the linear operator which appears is of third order. A good 
approximation can be obtained by discarding the term on 
the right hand side of eq.(29), which can then be readily 
integrated to give: 



BCr 

A 



rpr 



(38) 



In the above we have taken the integration constant which 
appears as being equal to zero, from requiring K- — > for 
Pr -> 0. Substituting the above relation for Vr into eq.(27) 
leads to: 



dpr 
dr 



a) 2 



(39) 



The second term in the above equation can be dismissed 
for r < 2GM/(?, which is in any case the validity regime 

defined by having previously neglected the right hand side 
of eq.(29), as it is easy to check from the final complete 
solution. In this regime we now obtain: 



Pr = Pr 



1/2 



/ \ 1/2 1/2 

Vr = — j Pr-r/ r I , 



(40) 



(41) 



Where rp is a characteristic radius at which pr = Pr- Now 
from eq.(28). 



(42) 



In the above equation we have also taken the integration 
constant as zero, from requiring Ur for r — s- 0. Choosing 
without loss of generality the two characteristic radii f and 
rp both equal to GMjc?, we can now write the full solution 
as: 



(34) Pi=p-.(^)"(3co.^(.)-l) 



^ (pA f2^V' ( icos\e) - 1 
' \p^)\Gm) \ gg 



\bpo )\Gm) \ 2^/^gg 



(43) 
(44) 
(45) 
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where we have introduced pj = pgpr and ge as the angular 
part of g{ff), g{0)/pQ. The full solution can be seen to de- 
pend only on the two parameters po and pj, normalisation 
constants for the densities of the background state and the 
perturbation, with the velocities of the perturbation solution 
depending only, and linearly, on the ratio Q = 2^^^pj/po, 
which is expected to be small, beyond natural dependences 
on the physical parameters of the problem, M and c. 

We see form eq.(45) that the angular velocity will be 
zero only for = 11/2, 9 = and 9 = n. Thus, movement 
along the plane of the disk will remain along the plane, but 
also, along the poles movement will be exclusively radial. 
This last point, together with the positive sign of the ra- 
dial velocity along the poles, provides for a well coUimatcd 
jet along the poles. From eq.(44) we sec that one has only 
to ask for a background state where matter is concentrated 
on the plane of the disk with relatively empty poles, e.g. 
g{9) for 6 — > and 9 ^ n, in order to obtain extremely 
large ejection velocities along the poles. The axial symmetry 
condition imposed on eq.(34) guarantees both axial symme- 
try and symmetry above and below the plane of the disk 
for the full solution. We see also that if one takes higher 
orders for m, the index of the Logondro polynomial solu- 
tion to eq.(33), one obtains increasingly more critical angles 
at positions intermediary between and 7r/2 at which the 
angular velocity goes to zero. In fact, more complex geome- 
tries and asymmetric jets appear, as inferred observationally 
by e.g. Fcrrcira et al. (2006), if m is taken as an arbitrary 
real number. However, modelling a situation where a polar 
jet dominates the ejection identifies m = 2 as the leading 
order. 

Notice that the qualitative behaviour of the solution 
is guaranteed by the exactness of the angular solution, the 
approximation r < 2GM/c^ used for solving the radial prob- 
lem will only introduce an error in the actual values of the 
radial velocities outside of r < 2GM/c^ , but will not change 
the fact that velocities will be of radial infall along the plane 
of the disk, 9 = 7r/2, and of radial outflow along the poles 
(where the background solution becomes very small), the jet 
solution for 9 = 0,it. 

The two constants of the problem, po and pj, can now 
be calculated once a choice of g{9) is specified, from the two 
conditions: 



Ma = 2n / poVosin{9)rd9, 
'0 



Mj = 2-K I piVisin{9)rd9, 



(46) 



(47) 



where 9j is a suitable angle defining the opening of the jet, 

in all likelihood very small, as will be see in the following 
section. In the above equations Ma gives the matter accre- 
tion rate onto the central star, and Mj the matter ejection 
rate due to the jet. Dimensionally, the two quantities above 
will scale as: 



Ma=Ca 



{nGMf _ 



Po 



po 



(48) 
(49) 



are two dimensionless constants which will 
depend on the choice of ge, and which would be expected to 
be of order 1. 

Qualitatively, this type of model naturally furnishes a 
tight disk-jet connection (c.f. cq. 44) e.g., as now firmly es- 
tablished in microquasars and AGN jets (see e.g. Marscher 
et al. 2002, Chatterjee et al. 2009). In the above systems 
bursts of enhanced jet activity are seen to follow temporal 
dips in disk luminosity output after small characteristic de- 
lay times. In the present model, such a situation would be 
expected if the critical radius for transition to radial flow in 
the disk made a sudden transition to higher values. Again, 
the drop in disk output might not reflect the disk material 
disappearing (in this case being swallowed by the central 
black hole, as sometimes proposed), but simply fading from 
view as heating mechanisms shut down, then naturally en- 
hancing jet activity as the effective Ma increases. 



4 PARTICULAR SOLUTIONS 

In order to present a sample of the trajectories expected 

in the model, we turn to the full solution to the problem, 
eqs.(21) and (22), but written in dimensionless form: 



dT 



d9_ ^ _^Q^i/2 ( sin{9)cos{9) 
dT 5 \ ge 



(50) 



(51) 



The above remain in spherical coordinates, where 71 = 
rc^/GM and T = t(? /GM. A choice of Q and ge now allows 
to numerically integrate trajectories. We take Q = 5 x 10~^, 
and 



ge 



(52) 



with 9o = 0.38 radians, or about 22 degrees. This values 
is far from representing an extremely flattened disk, and 
hence one which does not force the jet solution, see eq.(44). 
With those parameters, and initial conditions specifled in di- 
mensionless cylindrical coordinates as R = 2 and Z ranging 
from 0.4 to 1.4, we solve eqs.(50) and (51) through a flnite 
differences scheme to plot figure 1. For this case, the two 
lowermost curves present trajectories which all turn down- 
wards to converge onto the central star. These are solutions 
which essentially follow the background state, infalling onto 
the bottom of the potential well. As one raises the initial 
value of Z however, a threshold is crossed and curves of 
a very different type ensue, the six upper jet trajectories 
shown in figure 1. We see the large pressure gradients as- 
sociated with the distribution of matter in the background 
solution acting to break the fall of the incoming material, 
turn it back, and then accelerate it vertically through the 
vertical density gradients. These jet trajectories present a 
relatively constant thickening at the base of the jet of order 
R = 0.5, which then rapidly diminishes with height to even- 
tually yield an extremely collimated structure which rapidly 
narrows to below the resolution limit of the solution, notice 
the horizontal displacement of the vertical axis.. 

Note also that although the region beyond 7?. = 2 lies 
outside the approximation of the radial solution, the quali- 
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Figure 1. The figure gives trajectories for the complete solution, 
for a range of initial vertical height positions at i? = 2. The 
dimensionless parameter of the problem was chosen as Q = 5 X 
10~^. We see the lowermost two trajectories converging into the 
central potential, but the six upper ones shooting up into a polar 
jet, where velocities are over 4 orders of magnitude higher than 
what was originally present along the disk. 



tative form of the full solution will not deviate much from 
what is shown in figure 1, due to the exact character of the 
angular solution. The final well coUimated jet ejecting ma- 
terial along the poles will be a generic feature. 

We can also measure the velocities along the jet, which 
scale linearly with the value of Q, and are sensitive to the 
choice of gg and in this case Oq. For the cases presented 
here values upwards of lO^c commonly appear along the 
jet, much larger than the values of order c present in the 
disk. These radial velocities are somewhat inexact, but will 
always present large values for background state configura- 
tions which empty towards the poles (c.f. eq.44), and will 
present larger values as the vertical density gradients in- 
crease. Notice that the generic validity of the solution pro- 
posed does not require that all of the disk material should 
loose all of its angular momentum at the critical radius Rt, 
only that some of the disk material should loose most of 
its angular momentum at that point. Any minor, residual 
angular momentum remaining on the disk fraction which 
forms the jet will only establish a finite jet cross-sectional 
area through the appearance of a centrifugal barrier, which 
will present only a small correction on the scenario given 
here. Still, the empirical presence of spectroscopically stud- 
ied accretion disks truncated interior to certain critical radii 
suggests the very substantial reduction of the shears in that 
region, as will be the case when a transition to radial fiow 
takes place. 

In going back to eqs.(50) and (51) we see that one ex- 
pects 




Figure 2. The figure gives trajectories for the complete solution, 
for a range of initial vertical height positions at /? = 3.0. The 
dimensionless parameter of the problem was chosen as Q = 10~^. 
We see a very thin disk developing a central bulge inwards of 
R = 1, from which an extremely well coUimated jet of very fast 
material exits along the poles. 



Ma 



(53) 



It is reassuring that for the jets associated with T Tauri 
stars for example, values of Mj /Ma of between 0.1 and 0.01 
on average are observationally inferred and therefore of or- 
der 0.3 < Q < 0.1 (e.g Hartigan et al. 1996, GuUbring et 
al. 1998, Hartigan et al. 2004, Ferreira et al. 2006). These 
are compatible with the value of Q = 5 x 10"'^ used to plot 
fig. (1), and hence ones which will readily yield jet solutions. 
Given the substantial spread in the inferred values quoted 
above, we see that vales of 10""^ or even 10~* occur. We note 
that for smaller values of Q the jet becomes more centrally 
localised, making the numerical problem of integrating tra- 
jectories challenging. In general, thin disk jet solutions will 
appear for lower values of Q, or conversely, jet solutions for 
low values of Q require thinner disks, small values of ^o, for 
the function gg used in the example. 

Figure 2 presents an example with a much smaller value 
of Q = 10~*, requiring a slightly thinner disk with 9o — 0.24, 
close to 14 degrees. The figure shows a contour plot for di- 
mensionless momentum flow, the product of the dimension- 
less velocity, dlZ/dT, and density in units of M/Ti? from 
eqs. (43) and (44). Within the disk contours span the from 
—8 to —0.1, in intervals of 0.05, while in the jet contours 
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span from 2 x 10"* to 10~^ in intervals of 5 x 10"''. The 
disk appears clearly identified as being made up of infalling 
material, while the jet shows up as a region of pure outflow 
carrying little mass, but at velocities of order lO^c. Again, 
we sec the main features of the solution being well estab- 
lished in the region interior to 7i < 2, the details of the 
radial part of the solution beyond this region will be some- 
what off, within the qualitative solution shown. 

The long range stability and coherence of these struc- 
tures lies outside the scope of this work, and is in all prob- 
ability furnished by a series of mechanisms extensively ex- 
plored in the literature including angular momentum, mag- 
netic fields and pressure containment of the surrounding 
medium, e.g. Begelman et al. (1984), Blandford (1990), Falle 
(1991), Kaiser & Alexander (1997). 

In going to the more extreme jet phenomena associated 
with stellar black holes (e.g Mirabel & Rodriguez 1999), 
quasars (e.g. Marschor 2002) and gamma ray bursts (e.g. 
Mcszaros 2002), it is natural to expect the ideas presented 
here to apply, but amplified to a much more extreme regime 
by the appearance of corresponding relativistic and general 
relativistic effects, to first order, the shift in the divergence 
in the potential from r = in the Newtonian case to r = Vs 
in the general relativistic one. It is therefore natural to ex- 
pect purely hydrodynamic jet generation mechanism to ap- 
ply across all classes of astrophysical objects, specially given 
the qualitative scalings and similarities which appear over 
all astrophysical jet classes (e.g. Miarabel & Rodriguez 1999, 
Mendoza, Hernandez & Lee 2005), in addition to the mag- 
netically driven processes traditionally found in the litera- 
ture. 

Notice from eq.(44) the intimate link between the jet 
velocity and the physical state of the infalling material, c 
and Ma through po- This implies that temporal variations 
in the parameters of the infalling disk will result in temporal 
variations in the density and velocity of the jet material, in 
a way described by eqs.(43) and (44). The above can serve 
as a physical description of the key processes relevant to 
the formation of internal shocks in astrophysical jets, the 
main ingredient behind phenomena such as HH objects and 
gamma ray bursts. 



5 CONCLUSIONS 

We show that given a radially infalling accretion disk, a 
purely hydrodynamical jet ensues. 

We calculate the condition for the transition from quasi- 
circulax to quasi-radial flow in a standard accretion disk, 
and show it will always occur for power law surface density 
profiles of the form S oc interior to a critical radius, 

provided n > 1/2. 

Comparison with inferred inner holes in observed accre- 
tion disks yields results consistent with our estimates for the 
above transition radius, the point where shears in the fiow 
(and hence heating) end. 

Well collimated jets readily appear, proving the exis- 
tence of purely hydrodynamical mechanisms for the genera^ 
tion of astrophysical jets. 
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